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Abstract 

Background: Host-plasmid interactions liave been discussed largely in terms of the influences of plasmids, whereas 
the contributions of variations in host genomes to host interactions with foreign DNA remain unclear. A strain with 
a so-called "clean genome" {i.e., 1\/1DS42) of reduced genome size has recently been generated from the wild-type 
strain MGIGSS, a commonly used host strain. A quantitative evaluation of the influence of plasmid burdens in these 
two Escherichia coli strains can not only provide an understanding of how a reduced genome responds to foreign 
DNA but also offer insights into the proper application of these strains. 

Results: The decreases in growth caused by the cost of carrying foreign DNA were similar for the wild-type and 
clean-genome strains. A negative correlation between the growth rate and the total amount of exogenous DNA 
was observed in both strains, but a better theoretical fit with a higher statistical significance was found for the strain 
with the clean genome. Compared to the wild-type strain, the clean-genome strain exhibited a reduced carrying 
capacity for exogenous DNA, which was largely attributed to its ability to restrict the replication of foreign DNA. 
A tendency to allocate energy and resources toward gene expression, but not DNA replication, was observed in 
the strain with the clean genome. 

Conclusions: The possession of a clean genome constrained the plasmid copy number to a wild-type-equivalent 
load. The results indicate that the wild-type strain possesses a greater tolerance for foreign DNA, as in endosymbiosis, 
and that the use of strains with clean genomes will be favorable in the applications that require precise control and 
theoretical prediction. 
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Background 

Host-plasmid interactions have been intensively studied 
in both bacterial and yeast cells, with the aim of devel- 
oping a quantitative understanding of how exogenous 
DNA, i.e., plasmids, contribute to host growth [1-8]. De- 
creases in the growth rates of host cells that carry for- 
eign DNA {e.g., plasmids) have mainly been attributed to 
the additional costs associated with the expression of 
plasmid genes [9,10]. Because the size of the prokaryotic 
genome is considered to be constrained by bioenergetics, 
(gene expression across an entire genome is assumed to 
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consume approximately 75% of the total energy pro- 
duced by bacterial cells [11]), whether changes in the 
size of a host genome can improve a hosts capacity for 
DNA replication is an intriguing question. 

Furthermore, evaluations of the contributions of host 
genomes to host-plasmid interactions are also relevant 
to biotechnological applications. A strain with a clean 
genome of reduced size {i.e., MDS42) has recently been 
generated from the wild- type genome of strain MG1655 
[12]. The MDS42 genome possesses significant advan- 
tages in terms of plasmid stability [13,14], protein pro- 
duction [15] and other applications in synthetic biology 
[16,17]. These advantages for biotechnological applica- 
tions are largely attributed to genetic characteristics such 
as a low mutation rate and limited amount of recombin- 
ation [13,14,17]. 
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The latest genome-wide analyses have revealed that 
the reduced genome of MDS42 exhibits both a higher 
level of gene expression and a more fixed chromosomal 
periodicity [18] compared to its mother strain, MG1655 
[19]. These novel properties of MDS42 largely result 
from the complete deletion of insertion sequences [12], 
which are considered to be nonessential for cellular life. 
Thus, MDS42 is believed to be an ideal host for plasmid 
production [20] because plasmids exhibit a greater de- 
gree of structural stability in MDS42 compared to 
MG1655 [12]. However, to date, the increase in the 
productivity of plasmids inserted into MDS42 has not 
been well studied [21]. 

In this study, we investigated whether a reduction in 
host genome size would benefit the replication of foreign 
DNA by comparing two closely related strains {i.e., 
MG1655 and MDS42) of different genome sizes. Esti- 
mating the burden of foreign DNA {i.e., plasmids) in 
both MG1655 and MDS42 allows us to determine how a 
reduction in genome size influences the carrying cap- 
acity for foreign DNA, and this information also pro- 
vides further insight into the relationships among DNA 
replication, host growth and the cost distribution of cel- 
lular energetics. Any novel information concerning the 
nature of the host-plasmid interactions in strains with a 
clean genome {e.g., MDS42) can be taken into consider- 
ation when choosing strains for various applications. 

Results and discussion 

The growth cost of foreign DNA is similar for a wild-type 
or a clean genome 

To evaluate the contribution of the host genome to the car- 
rying capacity for foreign DNA, three plasmids of different 



sizes (S, M and L) were constructed, based on the se- 
quences of pUC19 and pSClOl (Figure lA). All three 
plasmids possessed a common replication site (on), a 
pMBl -derived origin from pUC19, and primarily followed 
identical replication mechanisms, in a temperature- 
inducible manner [22] . These plasmids were subsequently 
transformed into two closely related host strains of differ- 
ent genome sizes, the wild type E. coli strain MG1655 and 
its derivate MDS42, which possesses a reduced genome 
{i.e., a so-called clean genome'). Exponentially growing 
host cells bearing these plasmids were repeatedly analyzed 
at various temperatures (Additional file 1: Figure SI), to 
alter the plasmid copy numbers. Host cells that lacked 
plasmids were used as controls for data normalization at 
each temperature (Additional file 1: Figure S2). 

Decreases in growth rate were detected in both host 
strains when they contained plasmids (Figure IB, 
Additional file 1: Figure S2). A trend toward a decrease in 
host cell growth rate with an increase in plasmid size was 
observed. The decrease in growth rate caused by the pres- 
ence of foreign DNA was more significant when the host 
cells were grown at normal and higher than normal tem- 
peratures but was insignificant at the lowest temperature. 
The results strongly suggest that the burden of foreign 
DNA was roughly equivalent for the two hosts, although 
the genome of MDS42 was shorter in length. The costs as- 
sociated with the replication and/or expression of exogen- 
ous DNA were similar for the wild-type and the clean 
genomes. In addition, the results suggested not only the 
plasmid copy number (as previously reported [3,4]) but 
also the plasmid size (Figure IB) played a role in host 
growth. The total DNA amount, comprising the informa- 
tion of both copy number and plasmid size, might be a 
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Figure 1 Contributions of plasmids to cell growth. A. Schematic drawing of tine plasmids. All plasmids (pUC-S, pUC-M and pUC-L) were 
derived from pUC19. The lengths of S, M and L represent the plasmid sizes of 2,334, 6,426 and 9,573 bp, respectively. The site of initiation for 
replication and the ampicillin resistance gene are indicated as ori and bla, respectively. B. Cell growth rates. The exponential growth rates of E. coli 
strains MG1655 (upper panel) and MDS42 (lower panel) in a minimal medium are shown for different temperatures. The green, blue and red 
colors indicate the cells bearing the plasmids S, M and L, respectively. The strains that did not contain plasmids, which were used as controls 
{i.e., native growth fitness), are shown in black. The bars represent the standard errors of triplicates. Statistic significance is indicated (*, p < 0.05; 
p < 0.005). 
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proper parameter to evaluate the host-plasmid relations in 
common. 

The decrease in growth rate is dependent on the total 
amount of exogenous DNA 

In exponentially growing cells, the growth rates of the host 
cells depended more on the total number of nucleotides 
contained within the plasmids than on the number of 
plasmid molecules. Individual analyses of each plasmid- 
host construct revealed a highly significant negative cor- 
relation between normalized growth rates and plasmid 
copy numbers (PCNs) (Additional file 1: Figure S3). Nega- 
tive correlations were detected for all three plasmids and 
for both strains. This finding agrees well with the results 
of previous studies that examined the costs of plasmids to 
their hosts [9,10]. Intriguingly, even stronger correlations 
were found between host growth rates and the total length 
of exogenous DNA (Nt) when all three plasmid-host inter- 
actions were plotted together (Figure 2), indicating that 
the cost of plasmid replication for host growth is consider- 
able under the experimental conditions examined. The 
negative correlation between growth and PCN found 
when combining the data for all three plasmids was 



significant for MDS42 but not for MG1655 (Figure 2A). 
The higher significance and greater strength of the nega- 
tive correlation for M, compared to PCN, was a finding 
common to both genomes. Both the wild-type and the 
clean genomes appear to be similarly affected by the bur- 
den of foreign DNA. 

Enhanced predictability of the clean genome 

To achieve a quantitative understanding of the negative 
correlations described above, the adherence of the data 
to a simple linear regression model (model 1) and two 
models presented in previous studies [23,24] (models 2 
and 3) was assessed. The experimental data sets were 
used to fit these three models individually (Figure 3), 
and the parameter constants for these models were cal- 
culated by applying least squares regression (Table 1). 
The significances of the models were evaluated by their 
residual sum of squares (RSS) (Table 2). When consider- 
ing the relationship between host cell growth and M, we 
found that none of the three models provided a better fit 
to the data than the other models. However, model 3 
produced the lowest RSS value when the effect of PCN 
on growth rate was considered. Additionally, for all three 
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Figure 2 Relationship between cell growth and the amount of exogenous DNA. A. Relationship between cell growth and the plasmid copy 
number. B. Relationship between cell growth and the amount of exogenous DNA. The growth rate was normalized to a linear scale. The 
standard errors and correlation coefficients are given. The left and right panels show the results for strains MG1655 and MDS42, respectively. The 
color scheme is as described in Figure 1. 
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Figure 3 Comparison among three models. A. Relationships between cell growth rate and either plasmid copy number. B. Relationships 
between cell growth rate and the total amount of exogenous DNA.. The black and red open squares represent the data for the host strains 
MG1655 and MDS42, respectively. The bars represent the standard errors of triplicates. The solid, broken and chained lines denote the regression 
lines for the three models represented by Eqs. 1 (model 1), 2 (model 2) and 3 (model 3), respectively. 



models, the analyses of the effect of M produced lower 
RSS values compared to those of PCN, regardless of the 
host strain. These theoretical analyses support our new 
findings (Figure 2) and show that a better quantitative 
evaluation of the cost of plasmids on host growth is 
achieved by using the relationship between host growth 
and the total amount of DNA (M), as opposed to the 
plasmid copy number {PCN), 

Interestingly, the clean genome exhibited better fitting 
irrespective of the regression model employed. The data 
for the clean genome fit the models better than the na- 
tive host genome, and the RSS values were lower for 
MDS42 than MG1655 for both the PCN and Nt data 



Table 1 Estimated parameter constants in the models 



Host strain 


X 


a 


c 


n 


K 


MG1655 


PCN 


-1.7e-4 


3272 


0.25 


44 


MDS42 


PCN 


-5.6e-4 


973 


0.48 


106 


MG1655 


Nt 


-7.2e-5 


7340 


2.87 


1.3el0 


MDS42 


Nt 


-1.3e-4 


4065 


0.97 


5223 



The parameter constants {a, c, n and K) for the three models examined (Eqs. 1, 
2 and 3, details in the Methods section) were calculated using least squares 
regression. PCN and Nt are the copy number of plasmids and the total amount 
of plasmid DNA (the number of nucleotides), respectively. 



(Table 2). Both the increased significance of the correl- 
ation between host growth and the amount of exogen- 
ous DNA (Figure 3, Table 2) and the higher correlation 
coefficients (Figure 2) indicate the highly precise control 
of either cell growth or DNA replication in MDS42. 
More predictable levels of performance may occur in 
hosts with fewer nonessential sequences in its genome 
{i.e., MDS42). The results from this study were consist- 
ent with previous findings that have indicated a highly 
regular, predictable and well-organized behavior for 
MDS42 at the genome-wide expression level [18,25,26]. 



Table 2 Residual sums of squares (RSS) 



Host strains 


Models 


RSS_PCN 


RSS_A/t 


MG1655 


Eq. 1 


0.26 


0.09 


MG1655 


Eq.2 


0.26 


0.09 


MG1655 


Eq.3 


0.18 


0.08 


MDS42 


Eq.l 


0.11 


0.04 


MDS42 


Eq.2 


0.12 


0.04 


MDS42 


Eq.3 


0.09 


0.04 



The RSS values for the regression analyses were used to evaluate the 
three models. 
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Reduced carrying capacity of the clean genome 

Because the linear regression model (model 1) examines 
the relationship between cell growth rate and the 
amount of exogenous DNA, the slope {a) represents the 
magnitude of the growth cost of exogenous DNA to the 
host. A greater decrease in growth was clearly observed 
in MDS42 compared to MG1655 (Table 1), indicating 
that the difference in the carrying capacity for exogenous 
DNA depended on the host genome size. For example, 
when the linear regression was normalized by the gen- 
ome size (Additional file 1: Figure S5), in the case of a 
50% decrease in fitness, MG1655 could bear an amount 
of exogenous DNA that was 1.5 -fold greater than its 
genome size, whereas MDS42 could only bear an 
amount of exogenous DNAs that was equivalent to its 
genome size. The clean genome might save cellular en- 
ergy and/or resources for its own genome and restrict 
energy and/or resources spent on foreign DNA. 

A similar limit to the carrying capacity of foreign DNA 
was also detected based on cell size, as estimated by flow 
cytometry. The relative size of the host cells was positively 
correlated with the total amount of exogenous DNA 



(Figure 4A) and negatively correlated with the growth rate 
(Figure 4B), independent of the host strain. That is, the 
cell sizes of both host strains increased when they con- 
tained foreign DNA. However, MG1655 was able to 
provide a larger amount of cellular space than MDS42 
(Figure 4, Additional file 1: Figure S4), which could dir- 
ectly influence the carrying capacity of foreign DNA. The 
wild-type genome seemed to be more receptive of foreign 
DNA, while the clean genome appeared to provide more 
authoritarian control over the cellular conditions for for- 
eign DNA. 

Restricted replication allows for higher expression levels 
in the clean genome 

Further analysis demonstrated that the greater decrease 
in growth with increasing plasmid load {i.e., PCN and/or 
M) in MDS42 compared to MG1655 was due to a 
slower replication rate of the plasmids (Figure 5). The 
total amount of exogenous DNA is determined by the 
equilibrium between production and dilution. As both 
strains exhibited approximately equal growth rates under 
identical conditions (Figure IB), the dilution rates in the 
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Figure 4 Relationships between cell size and cell growth and between cell size and the amount of exogenous DNA. A. Relationship 
between cell size and the amount of exogenous DNA. B. Relationship between cell size and cell growth. Cell size is presented as the mean value 
of FSC signals measured using flow cytometry. The left and right panels indicate the strains MG1655 and MDS42, respectively. The color scheme 
is as described in Figure 1. 
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two strains were equivalent. Therefore, the accelerated 
decreases in the amount of exogenous DNA in MDS42 
(Table 1) were caused by decreases in the production 
(replication rates) of the plasmids. A reduced replication 
rate of exogenous DNA in MDS42 was observed inde- 
pendent of both the growth temperature (Figure 5, 
upper) and the sizes of plasmids (Figure 5, bottom). The 
reduction in host genome size decreased the replication 
capacity for foreign DNA. This result indicates the possi- 
bility that the removal of genes from MG1655 influ- 
enced the replication of foreign DNA in this strain. As 
the previously reported genes specifically related to the 
DNA and protein productivity [21] were not comprised 
in the deletions, there might be some novel functions 



interfering with replication or expression among those 
deleted genes. In another word, the regions that were re- 
moved most likely had a buffering effect, releasing the 
stress caused by the replication of exogenous DNA and 
helping to maintain cellular homeostasis. Further sys- 
tematic tests on different types of genomes along with 
transcriptome analysis are required to acquire a solid 
conclusion in mechanisms. 

In addition, the expression of the foreign DNA could 
affect the growth burden to a similar extent as replica- 
tion. To examine the degree to which the presence of a 
redundant gene decreases growth rates, the plasmid 
pUC-G (Figure 6A), which is approximately identical in 
size to pUC-S but carries a nonessential excess gene 
igfp), was inserted into the host cells. For both strains, 
the growth rates of the host cells that contained pUC-G 
were much slower than those of host cells that con- 
tained pUC-S (Figure 6B). This result verified that the 
gene expression of plasmids increased the growth bur- 
den beyond the cost paid for the replication of the 
foreign DNA alone. Importantly, MDS42 was more sen- 
sitive to this increase in gene expression; the growth rate 
of MDS42 was approximately one-third of the growth 
rate of MG1655 in this experiment. The high expression 
level of the foreign DNA caused the significant decline 
in the growth rate of MDS42. An approximately three- 
fold increase in the GFP expression level was detected in 
MDS42 (Figure 6C), which was consistent with the de- 
crease in growth rate. 

Taken together, the plasmid carrying capacity was lower 
in MDS42, independent of the growth temperature or the 
plasmid size (Figure 5), whereas gene expression was def- 
initely higher in MDS42 (Figure 6C), a finding in agree- 
ment with previous reports [15,18]. Compared to the 
wild-type genome, the clean genome assigned a higher 
priority to the gene expression of foreign DNA than to its 
replication. Such trade-off like phenomena of increased 
protein production with decreased plasmid carrying cap- 
acity in the reduced genome could be also detected when 
using the rich medium, despite of the highly accelerated 
cell growth (Additional file 1: Figure S6). This finding 
indicates that removing the redundant genes from a host 
genome can save energy and resources for their direction 
toward both host genome expression [18] and the ex- 
pression of foreign genes (Figure 6C), but not DNA 
replication. 

Conclusions 

In summary, the burden of foreign DNA was evaluated 
in two E. coli strains with different genome sizes. The 
growth rates of both host strains were more negatively 
correlated with the total amount of exogenous DNA 
than the plasmid copy numbers. The growth decrease, 
which was mediated by plasmid replication, was more 
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significant and more theoretically predictable in the 
strain with the clean genome (MDS42) than that with 
the wild-type genome (MG1655). The greater decrease 
in growth that occurred in MDS42 was caused by a 
slower replication rate of foreign DNA. The results indi- 
cate that the MG1655 strain can bear more exogenous 
DNA but experiences an equivalent loss in growth fit- 
ness, whereas MDS42 can restrict the amplification of 
foreign DNA to maintain its level of growth fitness. Plas- 
mids that are inhibited by constraints in the cellular 
conditions for replication might benefit from the highly 
induced level of gene expression found in MDS42. These 
results indicate that the wild- type strain MG1655 is 
more suitable for an endosymbiotic application and 
possesses an advantage in terms of adaptation and 



evolution, whereas the engineered, clean-genome strain 
MDS42 more closely resembles a type of living machinery 
(a controllable cell) and is more suitable for the precise 
control and production of synthetic materials. 

Methods 

Plasmid construction 

The lacZ region of pUC19 was removed using the In- 
Fusion HD Cloning Kit (Clontech) and the pUC19_del_ 
lacZ_Fw and pUC19_del_lacZ_Rv primers. Both colony 
PGR and blue-white selection were performed to verify the 
deletion. The resultant plasmid had a length of 2,334 bp 
and was named pUC-S. A DNA sequence (4,092 bp) of 
pSClOl was amplified with the pSC101_seq_Fw and 
pSC101_seq_Rv primers. Using the In- Fusion HD Cloning 
Kit, the resultant fragment was inserted into plasmid 
pUC19 in the same region where pUC-S was inserted, by 
using the pUC19_del_lacZ_linearize_Fw and pUC19_del_ 
lacZ_linearize_Rv primers. This increased the plasmid 
length of pUC-M to a total length of 6,426 bp. A longer 
fragment (6,986 bp) from pSClOl that included the 
tetracycline-resistance gene (Tc^ was amplified using the 
pSC101_seq_Fw2 and pSC101_seq_Rv2 primers and then 
inserted into pUC19 as described above to create the pUC- 
preL plasmid. Subsequently, a 198-bp fragment was re- 
moved from the translation initiation region for Tc^ in 
pUC-preL using the primers pSC101_del_TcR_Fw and 
pSC101_del_TcR_Rv, which formed a 9,573-bp-long plas- 
mid, pUC-L. Thus, all three plasmids that were produced 
had a common site of initiation for replication derived from 
pUC19, and carried a single gene for ampicillin resistance. 

A GFP sequence that includes the promoter Ptet was 
amplified from the pBRgalKGR plasmid (a pBR322 
derivative) [26,27] using the pBR _gfp_kanR_Fw and 
pBR _gfp_kanR_Rv primers. The lacZ region of pUC19 
was also removed using the In- Fusion HD Cloning Kit 
(Clontech) but with a different primer set (pUC19_del_ 
lacZ_linearize_Fw and pUC19_del_lacZ_linearize_Rv2), 
resulting in a 2,334-bp-long fragment, the same length 
as pUC-S. The two fragments were ligated using the 
In- Fusion HD Cloning Kit. The kan^ gene sequence 
(2,866 bp) was removed using the pUC_del_kanR_Fw 
and pUC_del_kanR_Rv primers, which finally created 
the pUC-G plasmid. EGOS™ Competent E, coli DH5a 
cells (Nippon Gene) were used for the cloning and amp- 
lification of the plasmids. All the primers used in this 
study are listed in Additional file 1: Table SI. 

Cell culture and FCM analysis 

Cells were cultured in minimal media as described previ- 
ously [26]. E, coli cells that carried the desired plasmids 
were grown in media that contained an additional 
50 (ig/mL of ampicillin. The cell concentrations were de- 
termined using a flow cytometer (FACSCanto™II, Becton, 
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Dickson and Company) equipped with a 488-nm argon 
laser and a 515-545-nm emission filter (GFP), as previ- 
ously reported [26]. The following PMT voltage settings 
were applied: forward scattering (FSC), 420; side scattering 
(SSC), 400; and FITC (GFP), 375. The flow cytometry 
datasets for FSC and GFP, which measure relative cell 
sizes and gene expression levels, respectively, were 
analyzed using custom-made R scripts, as described 
previously [28]. Cell growth was measured during the 
exponential growth phase; the initial and final cell concen- 
trations were approximately 1-2 x 10^ cells/mL and 0.5- 
2 X 10^ cells/mL, respectively, during this stage of growth. 
The growth rate {fi} was calculated according from the ini- 
tial and final cell concentrations and the length of time be- 
tween initial and final measurements (8 hours) [18]. The 
growth data was normalized (v) according to the following 
equations, as described previously [7,29] : 

fi = ln(Ct/Co)/t 



Evaluation of various models 

A simple linear regression model was applied to the data 
for both strains, based on the foUowing equation (Eq. 1, 
assigned to model 1): 



f{x) = ax-\- b 



(1) 



Here, a represents the cost bias of plasmid replication 
and b is the growth rate of the native host cells without 
plasmids. For example, when the host does not contain a 
plasmid {x = 0), the growth rate is equal to b {i.e., b = 1). 
In addition to the simple linear model (Eq. 1), the model 
primarily used in this study, two other models (Eqs. 2 
and 3) have also been previously used [23,24] to estimate 
the relationship between plasmid copy number {x) and 
host growth rate (y): 



K 



K + x"" 



(2) 
(3) 



Plasmid copy number and total amount of exogenous 
DNA 

Plasmid copy numbers were determined by qPCR using 
the StepOnePlus™ Real-Time PCR System (Applied Bio- 
systems) and SYBR Premix Ex Taq™ II (Tli RNaseH Plus) 
(TaKaRa). The primers (Additional file 1: Table SI) used 
to amplify the host genome {dxs_Fw and dxs_Rv) and the 
plasmid DNA {bla_Fw and bla_Rv) were identical to those 
used in a previous report [30]. Real-time qPCR reactions 
were performed using the following program: incubation 
at 95°C for 30 s, followed by 40 cycles of 95°C for 5 s and 
60°C for 60 s. After the final cycle, a melting curve analysis 
with a temperature gradient of 0.5°C/s from 60°C to 95°C 
was performed. Aliquots of purified genomic DNA and 
plasmid DNA were used as the standards for the genome- 
specific and plasmid-specific qPCR, respectively. The plas- 
mid DNA and genomic DNA were purified using the 
QIAprep Spin Miniprep Kit (QIAGEN) and the Wizard 
Genomic DNA Purification Kit (Pr omega), respectively, 
according to the manufacturers instructions. The abun- 
dances of both plasmid DNA and genomic DNA in the 
E. coli cells were determined using standard curves 
based on the bla and dxs genes, respectively. The plas- 
mid copy number (PCM) and the total amount of ex- 
ogenous DNA (Nt) were calculated according to the 
following equations: 

PCN = the abundance of plasmid DNA / 
the abundance of genomic DNA 



In all three models (equivalent to Eqs. 1, 2 and 3), y is 
the normalized growth rate and x indicates either the 
copy number of plasmids or the total amount of plasmid 
DNA. The parameter c represents the maximum pos- 
sible number of plasmids per ceU, and the exponent m 
was set to 0.5 according to the recommendation of a 
previous study [23]. 
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